Coupling of nitrogen-vacancy centers in diamond to a GaP waveguide 
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The optical coupling of guided modes in a GaP waveguide to nitrogen-vacancy (NV) centers in 
diamond is demonstrated. The electric field penetration into diamond and the loss of the guided 
mode are measured. The results indicate that the GaP-diamond system could be useful for realizing 
coupled micro cavity- NV devices for quantum information processing in diamond. 
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In recent years much progress has been made toward 
the realization of quantum information processing (QIP) 
based upon nitrogen-vacancy (NV) centers in diamond. 
Long 350 /is ground state electron spin coherence times 
have been observed [1] , full electron spin control has been 
achieved using optically-detected magnetic resonance [2 , 
and electron- nuclear qubit transfer, necessary for long 
quantum memory times, has been performed [3 . How- 
ever, such demonstrations so far have involved manip- 
ulation only of isolated NV centers. For realization of 
large-scale QIP [4] or for quantum repeater systems [5 it 
will be necessary to connect NV centers together through 
"flying" qubit s such as photons. To achieve this, optical 
structures in diamond such as microcavities and waveg- 
uides are needed to enable transfer of quantum infor- 
mation between the electron spin of the NV center and a 
photon. Optical structures in diamond could also be use- 
ful in enhancing photon collection efficiencies in recently 
proposed diamond-based magnetometers [6]. 

To date, two approaches have been pursued to couple 
NV centers to an optical device. The first is to couple 
NV centers in a diamond nanoparticle to a silica micro- 
cavity [7]. Unfortunately, nanoparticle NV centers ex- 
hibit spectral diffusion and large inhomogeneous broad- 
ening [8 making any practical, scalable device difficult to 
engineer. Groups have also worked on fabricating optical 
structures out of diamond [9l [T0|. Although this sec- 
ond approach is promising, the coupling of NV centers 
to a diamond optical device has not yet been demon- 
strated. In this work, we fabricate optical waveguides 
from high-index GaP (nGaP=3.3) thin films on top of 
diamond (nD=2.4) and demonstrate the coupling of the 
guided light to NV centers in the top 100 nm of the di- 
amond. We characterize the coupling strength of the 
guided mode to the surface NV centers as well as the 
scattering loss of the waveguides and calculate the ex- 
pected NV coupled GaP ring-cavity performance based 
on these results. 

The initial GaP wafer consisted of a commercially 
available (IQE) MBE-grown 250 nm GaP epilayer on an 
800 nm Alo.8Gao.2P sacrificial layer on a GaP substrate. 
The GaP epilayer was first thinned by ECR dry etch- 



ing to 120 nm to increase the penetration depth of the 
evanescent field once the waveguides have been trans- 
ferred to diamond. Optical ridge waveguides ranging in 
width from 1-8 /im with a final ridge height of 50 nm 
were defined using standard optical lithographic and dry 
etching techniques (Fig{l^,b). A 200 nm PECVD sili- 
con nitride film was deposited on some of the waveguide 
samples. The GaP epitaxial membrane was then under- 
cut and removed from the GaP substrate and transferred 
to diamond in a procedure adapted from Ref. ^J. In 
this procedure, a 7:100 IIF(49%):water solution at room 
temperature was used to etch the Alo.8Gao.2P sacrificial 
layer. The measured etch rate was 51 /im/hour along the 
(100) planes. 

The diamond samples were Sumitomo high-pressure 
high-temperature grown diamond with a typical nitrogen 
content of 30-100 ppm. The samples were implanted ei- 
ther with 2 MeV alpha particles with a dose of 10^^ cm~^ 
(HPHTl) , 200 keV Ga ions with a dose of 3 x 10^^ cm-^ 
(HPHT2), or 30 keV Ga ions with a dose range of 
10^^ — 10^^ cm~^ (HPHT3) in order to create vacancies. 
After implantation the samples were annealed at 925C 
in an H2/Ar forming gas during which vacancies com- 
bined with already-present nitrogen impurities to create 
NV centers [12 . 

The observation of NV photoluminescence (PL) from 
the implanted regions when exciting at 637 nm through 
the waveguide confirms the evanescent coupling of the 
guided electromagnetic field to NV centers near the sur- 
face. In Fig.[T]i 637 nm light from a continuous- wave tun- 
able diode laser is end coupled through a high-NA lens 
into an 8 /im-wide waveguide on HPHTl at room tem- 
perature. The waveguide crosses a dark diamond sector 
as well as several implantation regions. A long wave- 
length pass filter blocks scattered excitation light from 
the imaging CCD, and only light from the NV phonon 
sidebands with wavelength A > 653 nm is imaged. At 
cryogenic temperatures, different sectors in HPHT di- 
amond often exhibit different NV linewidths and peak 
intensities. In Fig. [l^ the waveguide coupling was per- 
formed at 10 K and a large contrast in intensity between 
the two diamond sectors is observed when the laser is 
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FIG. 1: (a) Waveguide structure before epitaxial liftoff, (b) SEM image of an etched 2 /xm-wide waveguide before liftoff. ECR 
etch conditions: Ar/BCls/Cb gases with 15/10/2 seem flow, 200 W microwave power, 45 W RF power, (c) Waveguide structure 
after liftoff and adhesion to diamond, (d) Optical excitation of NV centers through an 8 /xm waveguide. Implantation regions 
defined by the TEM grid and the dark sector boundary are clearly observed in the NV PL. (e) Sector boundary observed in 
the waveguide-excited PL in HPHT2. Spectra taken with 532 nm excitation from the top of the waveguide show the 637 nm 
NV" and 575 nm NV° ZPL hues in both sectors. 



slightly detuned from the narrower NV ZPL line in the 
lower sector. 

The two parameters critical to waveguide-NV coupled 
devices are the NV-waveguide coupling strength and the 
waveguide scattering loss. The coupling strength de- 
pends on the per-photon electric field strength of the 
guided mode at the NV site, and good coupling re- 
quires that the gap between the GaP membrane and the 
diamond surface is minimized. In order to determine 
the coupling strength the following experiment was per- 
formed. First a 250 nm GaP membrane was transferred 
to sample HPHT2. 532 nm laser light was end-coupled 
to the GaP membrane (Fig. |2^) and spectra contain- 
ing both the NV PL and the GaP LO phonon Raman 
signal were collected from the top with a lOx objective 
for TE and TM excitation (Fig. (2)3). To account for 
the relative brightness between the NV and GaP Raman 
signals, spectra were also obtained with top excitation 
from which a normalization factor was obtained that ac- 
counted for the GaP membrane thickness, the NV im- 
plantation depth (approximately 100 nm for 200 keV Ga 
ions), and the relative GaP Raman scattering efficiency 
between top and side excitation [13^ . The membrane was 
then thinned via an ECR dry etch and the optical mea- 
surements were repeated. As seen in Fig. |2]3, the ratio 
between the NV and the GaP Raman signals increased 
as the membrane was thinned and the evanescent field 
strength in the diamond increased. For comparison, fi- 
nite element simulations (COMSOL) were performed for 
a 1 fim wide waveguide to determine the relative electric 
field intensity in the waveguide and in the first 100 nm 
of the diamond for the lowest order TE and TM modes 
(Fig. [2]3 inset). The experimental and theoretical results 
for the dependence of the electric field strength ratio in 
the diamond and GaP as a function of membrane thick- 
ness are shown in Fig. Good agreement between the 
experimental data and the theoretical model can be ob- 
tained if a 4-5 nm air gap layer between the GaP and 



diamond is included in the model. This small gap causes 
only a modest reduction in the electric field strength for 
the TE mode. Currently we are unable to investigate 
the uniformity and composition of the diamond/GaP in- 
terface directly and can only say it is greater than the 
1 nm diamond surface roughness measured by AFM on 
a representative Sumitomo HPHT sample. 

The second parameter critical to GaP-NV diamond de- 
vices is the loss in the GaP at 637 nm. We first measured 
this loss for waveguides on the original GaP wafer (simi- 
lar to Fig.jl^). Fully etched waveguides with a waveguide 
height of 250 nm were fabricated with 6 lengths ranging 
from 0.1 to 3.2 mm. A 637 nm laser was coupled into 
the end of each waveguide through a 0.95 NA objective 
and the scattered output was collected from the top us- 
ing a lOx objective. No dependence on waveguide width 
within the measurement uncertainty was observed indi- 
cating sidewall roughness is not the dominant source of 
loss. The measured loss was 46 ±3 dB/cm (63 ±5 dB/cm) 
for the TE (TM) mode. This sample was next thinned 
to 160 nm and remeasured. The TE mode loss was un- 
changed, however the TM loss increased to 80 ±6 dB/cm. 
The TM mode has a higher electric field intensity at the 
top and bottom waveguide surface than the TE mode 
which could account for the increased TM loss with the 
waveguide thinning. 

The loss for GaP waveguides once transferred to di- 
amond was tested by exciting from the end of a ridged 
waveguide and measuring the intensity of the NV PL 
collected from the top of the waveguide in a single dia- 
mond sector. The longest decay lengths were obtained 
for waveguides which included a silicon nitride cladding 
layer on top of the GaP. A representative decay curve 
for a 1 fim wide waveguide is shown in Fig. [3| Loss for 
the TE (TM) mode in this waveguide was 72 ± 8 dB/cm 
(232 ± 13 dB/cm). The increased loss after waveguide 
transfer indicates that scattering loss at the horizon- 
tal surfaces is greater in the membrane/diamond system 
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FIG. 2: (a) A TE-polarized 532 nm laser is coupled into the overhanging region of a 250 nm GaP membrane on diamond 
and NV PL is imaged on a CCD camera, (b) Representative spectra collected from the top of the membrane for membrane 
thicknesses of 256 and 128 nm for TM excitation, (c) The ratio of the electric field strength in the top 100 nm of the diamond 
and in the GaP waveguide for various membrane thicknesses. Red (Blue) dots: TE (TM) experiment. Dotted lines: Red (Blue) 
TE (TM) calculation for a GaP waveguide on top of diamond. Solid lines: Red (Blue) TE (TM) calculation which includes a 
4.7 nm air gap between the GaP and diamond. Inset: Simulated TE/TM electric field intensity in a 1 /xm-wide, 120 nm thick 
waveguide (no air gap). 
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FIG. 3: (a) NV PL from a 1 micron waveguide excited at 637 
nm on sample HPHT3. (b) Decay of the PL in the 1 micron 
waveguide due to scattering loss. The measurement length is 
limited to the length of the diamond sector. 



than on the original GaP wafer. An increase in scattering 
loss is to be expected due to the larger refractive index 
change at the GaP /diamond interface, surface roughness 
on the diamond, and any damage/defects caused by the 
45 hour HF etch. Defects due to the HF etch can be ob- 
served under the optical microscope with defect densities 
increasing with etch time and HF concentration. Such 
defects could be eliminated if the GaP wafer was first 
bonded to the diamond and a substrate etching technique 
was used in place of the epitaxial lift-off technique. 

Previous work has shown that small-mode- volume pho- 
tonic bandgap cavities with cavity quality factors Q up 
to 1700 are attainable in a GaP membrane undercut on 
the original GaP substrate [14]. The losses we measured 
for GaP waveguides on diamond correspond to a max- 
imum possible cavity quality factor Q = 27rnGap/^(^ 
{a is the 1/e decay length) of 20,000 (6000) for the 
TE (TM) mode. This is already high enough that it 
should be possible to realize a large enhancement of the 
NV center-photon coupling efficiency in a GaP-based mi- 
crocavity structure. Finite element simulations (COM- 
SOL) [15j [16] show that a GaP ring cavity with a Q of 
20,000 (limited by a and not radiation loss) is enough to 



significantly enhance the NV emission into the ZPL line. 
The total spontaneous emission enhancement factor, in 
the weak-coupling limit and neglecting dispersion, is [17] 
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where g is the single-photon coupling strength. Hi is the 
cavity photon decay rate, 7totai = 13 MHz [18^ is the to- 
tal spontaneous emission rate including emission into the 
phonon sidebands, V = 18 (A/ncap)^ is the traveling- 
wave mode volume obtained from the numerical simula- 
tion, E{r^y) is the electric field at the NV site, ^(rViax) 
is the electric field maximum, and 7zpl = 0.35 MHz (de- 
rived from [T9^) is the spontaneous emission rate into 
the ZPL line. We find that an Fse > 1 is possible 
for a ring diameter of 2.5 /im coupled to an NV cen- 
ter 20 nm beneath the diamond surface. This total Fse 
factor corresponds to a 40-fold increase in emission into 
the zero phonon line. Such an increase in ZPL emission 
would be valuable for repeat-until-success schemes based 
on the DLCZ protocol [20 to create entanglement be- 
tween two remote NV centers through interference and 
detection ^|5j. 



SUPPLEMENTAL MATERIAL 
Ring-cavity geometry 

Fig. |4] shows the geometry in the numerical simulation 
used to determine the radiation-limited Q and mode vol- 
ume of a GaP on diamond ring cavity. The following 
parameters were used: GaP waveguide height = 120 nm, 
GaP waveguide width = 300 nm, GaP ring diameter = 
2.5 jam. For better mode confinement the diamond is 
etched down 120 nm to form a ridge. The NV is placed 
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FIG. 4: Simulation of the electric field \E\ in a GaP/ diamond 
ring cavity. 

20 nm beneath the GaP/diamond interface at the field 
maximum. 



Derivation of the spontaneous emission 
enhancement factor 

The expression given in Eq. 1 for the total on- 
resonance spontaneous emission enhancement factor in 
the weak-coupling limit is similar to the Purcell formula 
^SE = 3QX^ /Att'^V [21^ except for some additional fac- 
tors. 

For a single atomic transition z, if we neglect disper- 
sion in the dielectric, the ratio gf /ji given by Weisskopf- 
Wigner spontaneous emission theory is, 

^ 3 eD|^(rNv)pcos^^, ^^ 

where E{f) is the electric field of the cavity mode of in- 
terest, e(r) is the dielectric constant, tny is the location 
of the NV center, cji is the frequency of the atomic tran- 
sition, and Oi is the angle between E{f^y) and the dipole 
moment of the atomic transition. The factor enters 
through the photon density of states in bulk diamond, 
the factor eo = n'^ enters through the per-photon elec- 
tric field strength in bulk diamond, and the integral in 
the denominator enters through the per-photon electric 
field strength of the cavity mode. Combining this with a 
standard definition of mode volume, 

^ Jd're{r)\E{r)\^ 

e(w)|^(r-;nax)P ' ^ ^ 

and using e(f;nax) = n^^p, 7^ = 7zpl, Q = ^il^, and set- 
ting cos^^ = 1 (optimal dipole alignment) we obtain the 
result in Eq. 1. The factor 7zPL/7totai is important for 
NV centers since < 5% of the total spontaneous emission 
is into the zero-phonon line. If we omit this factor, we ob- 
tain instead the enhancement factor for the spontaneous 
emission rate through the zero-phonon line alone. 
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